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Abstract 

Partial sequences of the mitochondrial COI gene of forty eight European and two Turkish population 
samples of Myzus cerasi from different winter hosts {Prunus spp.) were subjected to phylogenetic analyses. 
The analysed M. cerasi samples emerged as paraphyletic relative to a Myzus horealis sample used as an out- 
group, and formed two major clades in neighbor joining, maximum parsimony, maximum likelihood and 
Bayesian inference trees, corresponding to subspecies living specifically on Prunus avium and P. cerasus. 
Multivariate discriminant analysis (method of canonical variates) was applied to find out if morphologi- 
cal variation of samples correlated with mitochondrial COI and host plant information. Mean scores on 
the first two canonical variables clustered samples fully in accordance with their COI haplotypes and 
host plants confirming the existence of two morphologically similar winter host - specific subspecies of 
M. cerasi in Europe. No single morphological character enabled satisfactory discrimination between apterous 
viviparous females of the two subspecies. A three-character linear discriminant function enabled 92.37% 
correct identification of apterous viviparous females of Af. cerasi cerasi (n=118) and 93.64% of M. cerasi 
pruniavium (n=l 10). A key for the morphological identification of the two subspecies is presented and 
their taxonomic status is discussed. 
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Introduction 

Black cherry aphid Myzus cerasi (Fabricius, 1775) is reported to be a serious pest of 
cherries all over the world (Barbagallo et al. 1997, Blackman and Eastop 2000, Hol- 
man 2009) and its morphology, life cycle, host specificity and potential harmfulness 
have therefore been the subject of intensive studies (Ross 1918, Wimshurst 1925, 
Pokrovskyj 1932, Vereshchagina 1966, Fomicheva 1967, Karczewska 1970, Rakaus- 
kas 1984, Gruppe 1990, Cichocka and Goszczynski 2004). Nevertheless, the species 
level classification of black cherry aphids has not been satisfactorily resolved. The black 
shiny aphids inhabiting cherry trees were originally described as a single species, Myzus 
cerasi (Fabricius, 1775), but European populations inhabiting sweet cherry {Prunus 
avium) were later separated as Myzus pruniavium Borner, 1926. Borner's species has 
been accepted by some (Borner 1943, 1952, Heinze 1961), but synonymised with 
cerasi by others (Miyazaki 1971, Eastop and Hille Ris Lambers 1976, Remaudiere and 
Remaudiere 1997), while others have treated it as a subspecies of cerasi (Shaposhnikov 
1972, Favret 2014). Differences in host specificity of the two taxa have been docu- 
mented in experimental transfer studies (Karczewska 1970, Dahl 1968, Vereshchagina 
1966), showing M. cerasi cerasi as heteroecious, alternating between cherries (both 
Prunus cerasus and P. avium) and herbaceous hosts {Galium, Euphrasia, Odontites, 
Veronica). M. cerasi pruniavium Borner, 1926 difi-ers from the nominative subspecies 
in having P. avium as the only winter host (it is incapable of living on P. cerasus), and 
also has a somewhat difi-erent summer host list. Enzyme electrophoretic studies indi- 
cated a reduced gene exchange between the subspecies (Gruppe 1988a), perhaps due 
to differences in the phloem sap of sour and sweet cherry causing divergent selection 
(Gruppe 1989). Morphological characters for discrimination between the two taxa 
have been suggested (Heinze 1961), but none of these enabled satisfactory separation 
between viviparous females of M. c. cerasi and M. c. pruniavium when applied to in- 
dependent aphid material (Dahl 1968, Gruppe 1988b). The taxonomic status of the 
black cherry aphid in the western Palaearctic thus remains unclear (Blackman and Eas- 
top 2000, 2006, Lampel and Meier 2007, Holman 2009, Osiadacz and Halaj 2010). 
They are possibly members of a complex of cryptic aphid species that includes M. cerasi 
umefoliae (Shinji, 1924), which overwinters on P. mume in Japan with Artemisia as its 
summer host (Miyazaki 1971). In northern Europe there is also another potentially 
distinct taxon that lives all year without host alternation on herbaceous plants (Dahl 
1968), and to which the name M. cerasi veronicae (Walker, 1848) may be applicable. 

A similar species complex is found in the mealy plum aphid, Hyalopterus spp. 
(Lozier et al. 2008). Separation of three Hyalopterus species was eventually justified 
by their distinctness at molecular level (Lozier et al. 2008), but they still remain dif- 
ficult to separate by their morphological characters (Basky and Szalay-Marszo 1987, 
Blackman and Eastop 1994, 2000, 2006), including those used in the most recent 
identification keys (Lozier et al. 2008, Rakauskas et al. 2013). A similar recent case is 
that of birch- and oak-inhabiting aphid species of the genus Stomaphis Walker, 1870 
(Depa et al. 2012). 
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The M. cerasi complex has not yet been the subject of detailed molecular study, 
although certain DNA sequences have been isolated (Foottit et al. 2008, Valenzuela 
et al. 2007, Clements et al. 2000, 2000a). Preliminary data on the partial sequences 
of the mitochondrial COI gene support the existence of subspecies of M. cerasi (Vo- 
ronova et al. 2011). The aim of this study was to clarify the taxonomic status of the 
host-alternating taxa in the M. cerasi complex by a combined study of partial sequences 
of the mitochondrial COI gene and morphological characters of the European samples 
collected from different species of cherries. 



Material and methods 

Material studied 

Fifty population samples of apterous viviparae of black cherry aphids from nine Eu- 
ropean countries and Turkey were collected in 2004-2013, mostly from P. cerasus, P. 
avium, but also from three other Prunus species (Table 1). Twenty samples were used 
for morphology - based stepwise discriminant analysis. The remaining 30 samples were 
used for subsequent evaluation of the derived discrimination functions. Samples of 
Myzus borealis OssiSinmhson, 1959 iLndMyzuspersicae {Sulzer, 1776) (GenBank Acces- 
sion No AB506741) were used as out-groups for the phylogenetic analyses. 



DNA extraction, PGR amplification and sequencing 

For molecular analysis, a single aphid individual from one sampled plant was consid- 
ered as a unique sample. Total genomic DNA was extracted from a single aphid using 
the DNeasy Blood & Tissue kit (Qiagen), which involved at least a 2 h digestion of 
tissue with proteinase K. Partial sequences of mitochondrial COI gene were PCR- 
amplified using earlier published primers (Turcinaviciene et al. 2006). PCR amplifi- 
cation was carried out in a thermal cycler (Eppendorf) in 50 [il volumes containing 
2 [il genomic DNA, 5 |.il of each primer (10 \iM), 5 |^1 of PCR- reaction buffer, 5 |^1 
of dNTP mix (2mM each), 4-8 ^1 of 25mM MgCl, and 1.25 U of AmpliTaq Gold 
360 polymerase (5U/|^1) and ddH^O to 50 [il. The cycling parameters were as follows: 
denaturizing at 95 °C for 10 min (1 cycle), denaturizing at 95 °C for 30", annealing 
at 49 °C for 30" and extension at 72 °C for 30" (37 cycles in total), and a final ex- 
tension for 5 min (1 cycle). PCR products were subjected to electrophoresis on 2% 
Top Vision agarose (Fermentas, Lithuania), stained with GelRed and sized against a 
MassRuler Low Range DNA ladder (Fermentas, Lithuania) under UV light. PCR 
products were purified and sequenced at Macrogen Europe (Amsterdam, the Nether- 
lands) and Institute of Biotechnology of the Vilnius University (Vilnius, Lithuania). 
The amplification primers were also used as sequencing primers. DNA sequences for 
each specimen were confirmed with both sense and anti-sense strands and aligned in 
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Table I . Aphid material used in the present study. Samples used for morphology-based discriminant 
analysis are given in bold. 



Place, date, collection number; (number of individual apterae per sample) 


GenBank Accession No 


Myzus cerasi on Prunus cerasus L. 


Jieznas, Prienai distr., Lithuania, 2012.05.30, 12-25; (8) 


KF754311 


Alytus, Lithuania, 2012.05.30, 12-30; (8) 


KF754310 


Daugai, Alytus disU., Lithuania, 2012.05.30, 12-32; (8) 


KF754303 


Skirgiskes, Vilnius distr., Lithuania, 2012.06.05, 12-37; (8) 


KF754304 


Eisiskes, Salcininkai distr., Lithuania, 2012.06.13, 12-43; (7) 


KF754329 


Labanoras, Svencionys distr., Lithuania, 2012.06.19, 12-70; (8) 


KF754306 


Moletai, Lithuania, 2012.06.19, 12-72; (8) 


KF754314 


Kraujaliai, Moletai distr., Lithuania, 2012.07.10, 12-120; (8) 


KF754305 


Zickai, Moletai distr., Lithuania, 2012.07.13, 12-132; (8) 


KF754321 


Nida, Neringa, Lithuania, 2012.08.09, 12-176; (8) 


KF754322 


Smiltyne, Klaipeda, Lithuania, 2012.08.12, 12-191; (8) 


KF754309 


Preila, Neringa, Lithuania, 2012.08.13, 12-203; (8) 


KF754330 


Bagnolo Mella, Brescia prov, Italy 2013.05.01, 13-27; (8) 


KF754338 


Poncarale, Brescia prov., Italy 2013.05.02, 13-33; (8) 


KF754337 


Suginciai, Moletai distr., Lithuania, 2013.06.15, 13-83; (8) 


KF754345 


Akmeniai, Lazdijai distr., Lithuania, 2013.05.30, 13-57; (8) 


KF754341 


Karsava, Ludza distr., Latvia, 2013.07.17, 13-133; (8) 


KF754344 


Gorodok, Vitebsk distr., Belarus, 2008.06.17, 08-6; (6) 


KF754302 


Zadrachje, Vitebsk distr., Belarus, 2008.06.18, 08-18; (8) 


KF754325 


Riga, Latvia, 2008.07.03, 08-73; (8) 


KF754328 


Skirgiskes, Vilnius distr., Lithuania, 2011.06.15, 11-46; (8) 


KF754316 


Cluj Gilau, Romania, 2012.06.19, Z12-112; (8) 


KF754327 


Cluj Gilau, Romania, 2012.06.19, Z12-116; (8) 


KF754319 


Poncarale, Brescia prov., Italy 2013.05.02, 13-30; (8) 


KF754346 


Mezopeterd, Haidu-Bihar distr., Hungary 2012.06.20, Z12-122; (8) 


KF754333 


Myzus cerasi on Prunus avium L. 


Skirgiskes, Vilnius distr., Lithuania, 2012.06.05, 12-39; (7) 


KF754332 


Salcininkai, Lithuania, 2012.06.13, 12-48; (6) 


KF754317 


Bratoniskes, Vilnius distr., Lithuania, 2012.06.14, 12-56; (7) 


KF754313 


Blagoevgrad, Bulgaria, 2012.06.26, 12-83; (7) 


KF754318 


Frankfurt/Main, Germany 2012.06.30, 12-104; (8) 


KF754307 


Kraujaliai, Moletai distr., Lithuania, 2012.07.10, 12-111; (8) 


KF754308 


Stirniai, Moletai distr., Lithuania, 2012.07.12, 12-128; (8) 


KF754320 


Juodkrante, Neringa, Lithuania, 2012.08.10, 12-182; (8) 


KF754323 


Pervalka, Neringa, Lithuania, 2012.08.11, 12-188; (7) 


KF754336 


Preila, Neringa, Lithuania, 2012.08.13, 12-199; (7) 


KF754335 


Rondo, Katowice, Poland, 2011.05.13, 11-10; (8) 


KF754349 


Tekir, Karamanmarash distr., Turkey 2011.05.21, 11-25; (5) 


KF754312 


Goksun, Karamanmarash distr., Turkey, 2011.05.21, 11-27; (8) 


KF754315 


Kairenai, Vilnius, Lithuania, 2010.07.01, 10-3; (7) 


KF754324 


Zafferana, Catania, Italy 2004.06.28, 04-49; (5) 


KF754339 


Costinesti, Romania, 2012.06.15, Z12-90; (8) 


KF754326 


Galata, Varna, Bulgaria, 2012.06.18, Z12-102; (7) 


KF754351 


Burgas, Bulgaria, 2012.06.19, Z12-110; (8) 


KF754331 


Cluj Gilau, Romania, 2012.06.19, Z12-117; (7) 


ICF747679 
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Place, date, collection number; (number of individual apterae per sample) 


GenBank Accession No 


Carpendolo, Brescia prov., Italy, 2013.04.27, 13-12; (7) 


KF754340 


Akmeniai, Lazdijai distr., Lithuania, 2013.05.30, 13-60; (8) 


KF754347 


Wojslawice, Lower Silesia, Poland, 2013.06.20, 13-98; (8) 


KF754342 


Myzus cerasi on Prunus serrulata Lindl. 


Wojslawice, Lower Silesia, Poland, 2013.06.20, 13-97; (8) 


KF754348 


Myzus cerasi on Prunus maackii Rupr. 


Dobele, Latvia, 2013.07.03, 13-119; (8) 


KF754343 


Myzus cerasi on Prunus mahaleh L. 


Medias, Sibiu distr., Romania, 2012.06.19, Z12-1 13; (8) 


KF754334 


Myzus borealis on Galium rubioides L. 


Zmejinyje ostrova, Kanev distr., Cherkasy reg., Ukraine, 2006.06.16, 06-74 


KF754350 



the BioEdit Sequence Alignment Editor (Hall 1999). Partial sequences were tested 
for stop codons and none were found. The sequence data have been submitted to the 
GenBank, accession numbers are given in Table 1 . 



Analysis of DNA sequences 

In addition to the sequences from 50 samples oiM. cerasi, COI sequences of Af. borealis 
from subgenus Myzus sensu stricto (the same subgenus as M. cerasi) and M. persicae 
from subgenus Nectarosiphon Schouteden, 1901 were selected as out-groups for the 
phylogenetic analyses, which included neighbor joining (NJ), maximum parsimony 
(MP), maximum Ukelihood (ML) and Bayesian inference in phylogeny (BI). NJ, MP 
and ML analyses were performed using MEGA 5 (Tamura et al. 201 1). For NJ analysis 
Kimura 2-parameter (K2P) model of base substitution was used. Bootstrap values for 
NJ, MP and ML trees were generated from 1000 replicates. For ML analysis Tamura 
3-parameter model with Gamma distribution (T92+G) was selected by MEGA 5 mod- 
el selection option (Tamura et al. 201 1). Bayesian analysis was conducted in MrBayes 
3.2.1 (Ronquist and Huelsenbeck 2003) using Hasegawa-Kishino-Yano model with 
Invariable sites and Gamma distribution (HKY+I+G), which was selected by jModeltest 
(Posada 2008). Four simultaneous chains, 3 heated and 1 "cold", were run for 3 000 
000 generations with tree sampling every 1000 generations. The topologies obtained by 
NJ, MP, ML and BI were similar, so only ML tree is shown with values of NJ/MP and 
ML/BI bootstrap support and posterior probabilities over 50% indicated above and 
below branches respectively Statistical parsimony haplotype networks were constructed 
for samples oiM. cerasi and Af. borealis using TCS v 1.21 (Clement et al. 2000). 



Morphological study and discrimination analysis 

Samples representing different clades in the molecular tree and haplotype network 
were used for stepwise discriminant analysis followed by canonical analysis: 10 sam- 
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pies from sour cherry, P. cerasus, and 10 samples from sweet cherry, P. avium (shown 
in bold in Table 1). 

Based on earlier taxonomic work (Heinze 1961, Dahl 1968), 19 metric (in mm) 
characters were studied: A2L - length of antennal segment 2; A3L - length of antennal 
segment 3; A4L — length of antennal segment 4; A5L — length of antennal segment 
5; A6BL — length of basal part of antennal segment 6; A6TPL — length of terminal 
process of antennal segment 6; BL - body length (excluding cauda); Bwant3 - basal 
width of antennal segment 3; CL - length of cauda; CW - basal width of cauda; DT3L 
- length of the second segment of hind tarsus; F3L - length of hind femur; FF - depth 
of the frontal furrow; SL - length of siphunculus; T3L - length of hind tibia; URL - 
length of ultimate rostral segment; URW - basal width of ultimate rostral segment; 
VBSLmax — maximal length of the ventral body hairs; VBSLmin — minimal length of 
the ventral body hairs. 

Measurements of the slide-mounted apterous viviparous females were performed 
by means of interactive measurement system Micro-Image (Olympus Optical Co. 
GmbH). STATISTICA 8 version software (Statsoft 2007) was exploited for data 
analysis. Pearson's correlation coefficients were calculated to evaluate the correlation 
of morphometric characters with body length. Characters with strong (| r | > 0.70) 
statistically significant (p<0.05) correlation with body length were removed from the 
further analysis: BL (r=1.00), F3L (r=0.83), T3L (r=0.82), A2L (r=0.7), A3L (r=0.75), 
A4L (r=0.71), A5L (r=0.7). The remaining 12 characters were used for the forward 
stepwise discriminant analysis followed by canonical analysis, with sample collection 
number as the grouping variable, thus excluding information about the host plant 
from the analysis. Mean canonical scores of the first two canonical variables were rep- 
resented as bivariate scatter plots, in order to show any clustering of samples. 

Morphological characters that contributed most to canonical discrimination func- 
tions were evaluated as having potential for separation of taxa. An identification key was 
constructed based on these discrimination functions. The key was then tested on the 30 
aphid samples that had not been used in its construction (listed in normal font in Table 1). 



Results 

Partial sequences of mitochondrial COI gene 

Fifty partial COI sequences of M. cerasi and one of M. borealis from 1 1 countries were 
included in analysis. The alignment contained 616 bases in the final set with three 
variable sites, all of which appeared parsimony informative. The average base composi- 
tion was A = 34.0%, C = 12.7%, G = 12.3% and T = 41.0%. The overall transition/ 
transversion ratio R = 1.221 for all sites. 

Five COI haplotypes were detected (Fig. 1): one for M. borealis, two for samples 
from P. cerasus and two for samples from P. avium (Table 2). Aphids collected from 
P. mahaleb, P. maackii and P. serrulata had the same haplotype (No. 3) as the majority 
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Figure I. Haplotype nerwork (TCS 1.21 software: Clement et al. 2000) for COI fragment (616 posi- 
tions in final set) haplotypes of Myzus cerasi and Myzus boreaiis. The haplotype with the highest outgroup 
probability is displayed as a square, while others are displayed as ovals. For sample information, see Table 2. 



Table 2. COI haplotypes of three Myzus taxa revealed by construction of haplotype network. Sample 
numbers are the same as given in Table 1 . 



Haplotype 
number 


Number of 
sequences 


Sequence 
length (bp) 


Sample numbers 


M. cerasi cerasi (collected from Prunus cerasus) 


1 


2 


616 


13-33; 13-27. 


2 


23 


616 


08-6;12-32; 12-176; 12-30; 08-18; zl2-122; 12-120; zl2-l 12; 12-132; 
12-25; 12-43; 12-191; 12-203; 12-70; 12-72; 08-73; 13-83; 13-133; 11- 
46; zl2-l 16; 12-37; 13-30; 13-57 


M. cerasi pruniavium (collected from Prunus avium, except where otherwise noted) 


3 


23 


616 


11-10; 12-39; zl2-l 10; 12-182; 12-104; 12-56; 12-83; 12-111; 12-199; 
zl2-90; Z12-102; 12-48; 12-188; zl2-l 13 (P. maackii); 04-49; 12-128; 
10-03; 13-60; 13-98; 13-97 (P. serrulata); 13-12; 13-119 (P. mahaleh); 
Z12-117 


4 


2 


616 


11-27; 11-25 


M. boreaiis (collected from Galium rubioides) 


5 


1 


616 


06-74 



of samples from P. avium. COI haplotypes detected among samples from P. cerasus 
(No. 1 and 2) and the remaining Prunus species (No. 3 and 4) differed in the follow- 
ing nucleotide positions of 616 bp ahgnment: 300 (between No. 1-2 and No. 3-4), 
321 (between No. 3 and No. 4) and 390 (between No. 1 and No. 2). The range of the 
intraspecific pairwise sample divergences (K2P model) was 0.0 - 0.5% (average 0.2%), 
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Table 3. Range of pairwise interspecific sample divergences of mitochondrial COI gene fragment (K2P 
model) for three species oi Myzus (number of samples used is in parentheses). 



Species 1 


Species 2 


Range of divergence, % 


M. cerasi (50) 


M. borealis ( 1 ) 


0.2-0.5 


M. cerasi (50) 


M. persicae ( 1 ) 


6.6-6.8 


M. borealis (1) 


M. persicae ( 1 ) 


6.8 



whilst interspecific pairwise sample divergences between three species of Myzus ranged 
from 0.2 to 6.8% (Table 3). 

The maximum parsimony (MP) analysis of partial COI sequences resulted in 930 
equally parsimonious trees (length = 43, CI=1.00, RI=1.00). The ML tree (T92 model) 
showed similar topology, as did NJ (K2P distances) and BI (HKY+I+G model) analyses. 
NJ, MP and ML bootstrap values over 50% together with BI posterior probabilities 
over 0. 50 are given at respective nodes of the same tree in Fig. 2. Thus the M. cerasi sam- 
ples form two major clades corresponding to two host-specific black cherry aphid taxa. 
One clade consists of all but two of the samples from P. avium, plus aphids collected 
from P. mahaleb, P. maackii and P. serrulata. The other clade contains all samples from 
P. cerasus and also includes the sample of M. borealis collected from Galium rubioides. 



Morphology 

When morphometric data of apterous viviparous females from 20 difi-erent geographi- 
cal localities were subjected to discriminant analysis with sample collection number 
as the grouping variable, the first two canonical variates (Fig. 3) clearly separated sour 
cherry samples (COI haplotype No. 2) from those collected from sweet cherry (COI 
haplotype No. 3). Length of terminal process of antennal segment 6 (A6TPL), length 
of siphunculus (SL) and maximal length of the ventral body hairs (VBSLmax) ap- 
peared to be important predictors for separation of the two taxa (Table 4). 

To discriminate between apterous viviparous females of host-specific black cherry 
aphid samples representing different clades in the haplotype network and the phyloge- 
netic tree (Figs 1-2), the following linear discriminant function (LDF) was obtained: 
3.924682xSL - 5.6667xA6TPL - 32.5504xVBSLmax + 1. Using this LDF, 97.37% 
individuals from the whole dataset were reclassified correctly into their a priori speci- 
fied groups with host plant species as grouping variable, including 96.2% of apterous 
viviparous females from P. avium (n=79) and 98.6% from P. cerasus (n=73). The post 
hoc classification of the remaining thirty samples gave 92.37% correct specimen iden- 
tification of M. cerasi cerasi (n= 118) and 93.64% of M. cerasi pruniavium (n= 110). The 
scatterplot of the mean LDF and body length values calculated for each of 30 samples 
representing difi-erent host specific subspecies of M. cerasi is shown in Fig. 4. The fol- 
lowing key is therefore suggested for the identification of apterous viviparous females 
of the two subspecies oi Myzus cerasi when sampled from winter hosts. 
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- P. cerasus LT 12-120 

- P. cerasus LT 12-132 

- 'P. cerasus LT 12-1 76 

- P. cerasus BY 08-1 8 

- *P. cerasus ROz12-112 

- P. cerasus LV 08-73 

- * P. cerasus LT 1 2-43 

- *P. cerasus LT 12-203 

- P. cerasus ROz12-116 

- P. cerasus HUz12-122 

- P. cerasus LT 13-57 

- P. cerasus LV 13-133 

- P. cerasus LT 13-83 

- P. cerasus IT 13-33 

- P. cerasus IT 13-27 

- P. awum TR 11-25 

- P. aWum TR 11-27 

-Myzus [Myzus) boreaiis Galium rubioides UA 06-74 
-IVIyzus {Nectarosiphon) persicae AB506741 



Figure 2. Maximum likelihood (ML) tree showing phylogenetic relationships among Myzus cerasi based 
on partial sequences of mitochondrial COI (616 positions in final set). Numbers above branches indicate 
support of NJ (left, > 50%) and MP (right, > 50%) bootstrap test with 1000 replicates, and numbers below 
branches indicate support of ML (left, > 50%) bootstrap test with 1000 replicates and posterior probabili- 
ties of BI analysis (right, > 0.50). Samples used for the discriminant analysis with a priori specified group 
membership followed by the construction of identification key are asterisked (*) . The remaining samples 
were used for the post hoc classification. Sample numbers are the same as given in Table 1 , together with 
the abbreviated symbol of respective country BG - Bulgaria, BY - Belarus, D - Germany, HU - Hungary, 
IT - Italy, LV - Latvia, LT - Lithuania, PL - Poland, RO - Romania, TR - Turkey, UA - Ukraine. 
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Figure 3. Plot of the mean scores of the first two canonical variates for 20 samples of Myzus cerasi (for 
specimen numbers per sample see Table 1). Samples cluster in accordance with winter host plant and COI 
haplotype (haplotype number is given in parentheses, see Table 2 for other haplotypes) . 
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Figure 4. Plot of the mean scores of the individual LDF values (number of specimens per sample is given 
in Table 1) plotted against the mean body length for 30 samples of Myzus cemsi (normal font in Table 1) 
used to evaluate eflictiveness of the eventual identification key. The icons are color-coded to match the 
COI haplotypes. Samples cluster in accordance with winter host plant and COI haplotype (haplotype 
number is given in parentheses, see Table 2 for haplotype information). 
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Table 4. Contributions of 1 1 morphological characters to the canonical function discriminating 20 sam- 
ples of M. cerasi. Character abbreviations are the same as in the text (Material and methods). 





WilKS i_,ainDCLa 
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U.UU 
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Ad O/^ 


A AA 


U.31 
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U.2 i 


Dwantj 




u.v 1 




A AA 
U.UU 


n 71 
U./ i 


A 90 

u.zv 


URL 


0.23 


0.93 


9.51 


0.00 


0.65 


0.35 


A6BL 


0.23 


0.94 


8.69 


0.00 


0.47 


0.53 


DT3L 


0.23 


0.97 


4.62 


0.01 


0.49 


0.51 


FF 


0.23 


0.97 


4.23 


0.02 


0.86 


0.14 


VBSLmin 


0.23 


0.97 


3.70 


0.03 


0.85 


0.15 


CL 


0.22 


0.98 


2.38 


0.10 


0.57 


0.43 


CW 


0.22 


0.99 


1.75 


0.18 


0.56 


0.44 



Key to European subspecies of Myzus cerasi on Prunus (apterous viviparous females) 

1 Value of LDF [3.92x( length of siphunculus) - 5. 67x (length of terminal process 
of antennal segment 6) - 32.55x(maximal length of the ventral body hairs) +1] 
greater than zero. On P. cerasus (and sometimes on P. avium) M. c. cerasi 

— Value of LDF less than zero. On P. avium, P. maackii, P. mahaleb, P. serrulata... 

M. c. pruniavium 



Discussion and conclusions 

The combination of genetic distance evaluation with phylogenetic tree-building meth- 
ods and multivariate analyses of morphometric data has been successfully applied to 
solve taxonomic problems in aphids, particularly in the genera Hyalopterus (Lozier et al. 
2008), Pentalonia (Foottit et al. 2010), Aulacorthum and Neoaulacorthum {=Pseudomeg- 
oura) (Lee et al. 201 la). Based on the global data set, the average genetic divergence of 
COI barcode sequences between aphid species within the same genus was reported to 
be 5.84% (range 0.46 - 11.3%), and that within species 0.05% (0.00-1.00%) (Foot- 
tit et al. 2008, Lee et al. 2011b). Interspecific divergence of six species representing 
three subgenera of Myzus calculated for COI barcode sequences was reported as ranging 
from 5.55 to 11.3% (Foottit et al. 2008). In comparison, partial COI sequences (Gen- 
Bank, 1145 bp) of three yl/)fe^^<2d" subspecies {A.f.fabae, GenBank accession numbers 
FJ965713, FJ965717-FJ965718; A/ cirsiiacanthoidis, FJ965698-FJ965709; A/ mor- 
dvilkoi, FJ965710-FJ965712) show values of genetic divergence (K2P model) ranging 
from 0.00 to 1.42% (Table 5). Therefore, the range of genetic divergence between the 
two clades of M. cerasi emerging in phylogenetic trees presented in this paper (0.0 to 
0.5%) appears to be of intraspecific level. Based on the available COI data, black cherry 
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Table 5. Pairwise sample divergences of 1145 bp mitochondrial COl gene fragment (K2P model) be- 
tween three subspecies o{ Aphis fabae (number of sequences used is in parentheses). 



Subspecies 1 


Subspecies 2 


Mean and range of divergence, % 


AJ.fabaelS} 


A. f. cirsiiacanthoidis (12) 


1.2 (0.97-1.42) 


A. f. fabae (i) 


A. f. mordvilkoi (3) 


0.15 (0.00-0.26) 


A.f. cirsiiacanthoidis (12) 


A. f. mordvilkoi (3) 


1.05 (0.97-1.15) 



aphids inhabiting sour and sweet cherries should therefore still be regarded as a single 
species. 

M. borealis is clearly closely related to M. cerasi and differs by only 0.2-0.5% of 
the COI sequences involved in the analysis. This suggests that it may also belong to 
the same species level taxon. More samples of M. borealis are needed to confirm this 
hypothesis. However, it should be noted that partial COI sequences of two biologically 
distinct Macrosiphum species, M. rosae (Linnaeus, 1758) and M. knautiae Holman, 
1972 are very similar (Turcinaviciene and Rakauskas 2009), and low divergence levels 
have also been reported for Bursaphis species (Rakauskas et al. 2011) and adelgids 
(Zurovcova et al. 2010). It seems probable that in rapidly speciating aphid groups one 
may expect to find low levels of COI sequence divergence between taxa that are never- 
theless functioning as distinct species. 
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